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Q I In the framework of the unimodular metagravity, with the scalar graviton/graviscalar dark 

matter, a regular anomalous one-parameter solution to the static spherically symmetric 
metagravity equations in empty space is found. The solution presents a smooth graviscalar 
halo, with a finite central density profile, qualitatively reproducing the asymptotically fiat 
CJ ! rotation curves of galaxies. To refine the description studying the axisymmetric case in 

the presence of luminous matter is in order. 
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1- Introduction 

O I General Relativity (GR) is well-known to be a theory of the massless tensor graviton, 

with two physical degrees of freedom residing in the metric field. At that, the unphysical 
metric components are restrained due to the general covariance (GC). In this respect, the 
unimodular covariance (UC) is a viable alternative to GC. Namely, it has been shown that 
rS I UC is necessary and sufficient to retain two transverse components for a massless tensor 

I field, with GC being thus excessive to this purpose.^ This may de facto justify all the GC 

violating alternatives to GR, possessing the residual UC. Viz., they may be considered 
as the theories of the massless tensor graviton with the different realizations of a metric 
component corresponding to a dilaton. Retaining this component but making it unphysi- 
cal due to GC we would arrive at GR. Respectively, two routs to go beyond GR, with the 
residual UC, are envisaged, (i) To eliminate such a component from the metric ab initio 
by means of a unimodularity condition. This would bar the local scale transformations, 
with the local measure becoming an absolute element.^ This is the so-called unimodular 
relativity /gravity. (For a recent discussion, see, e.g., Refs. 3 - 4.) The cosmological con- 
stant emerges here as an integration constant, instead of a Lagrangian parameter in GR. 
Thereupon, one hopes to naturally explain the cosmological constant being tiny (the long- 
standing naturalness problem). Furthermore, extending the unimodular relativity/gravity 
by an exactly massless dilaton one can try to explain the hierarchy problem in SM unified 
with gravity and simultaneously solve the so-called dark energy problem.^ (ii) To convert 
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the aforesaid unphysical (but still "harmless") component to the physical one by adding 
to the GR Lagrangian a GC violating term with a derivative of the metric. At that, 
a (massive) dilaton would part of the metric field.® Refiecting GC violation, 

such a rout originally implies (a class of) the distinguished, "canonical" coordinates, with 
the restricted (unimodular) group of the admitted transformations (thereof, the so-called 
restricted relativity/gravity). 

In an earlier paper, we put forward a hypothesis that GC violation with the metric 
derivative terms may serve as a raison d'etre for appearance in the Universe of the dark 
matter (DM) of gravitational origin.^ The reason is that under such a GC violation the 
energy-momentum tensor of the ordinary matter alone ceases to be covariantly conserved. 
This non-conservation can be compensated by equivalently treating the additional terms 
in the gravity equations as an energy-momentum tensor for the additional gravity degrees 
of freedom (the extra "gravitons") and associating the latter ones with DM. The metric 
itself serves thus as a resource of DM. As a simplest realization of this approach, the 
residual UC, with the local scale covariance alone being violated, was imposed. In this 
case, the metric comprises just one extra physical degree of freedom, a (massive) scalar 
graviton/graviscalar besides the (massless) tensor graviton. Possessing UC and containing 
the extra graviton, such a theory may be called the "unimodular metagravity" . By 
introducing a non-dynamical scalar density we put the theory to the arbitrary observer's 
coordinates, beyond the canonical ones. In a subsequent paper, the graviscalar field was 
taken as an independent variable substituting a metric component in the desired observer's 
coordinates.^ This allowed us to straightforwardly confront the unimodular metagravity 
with GR in the presence of an ordinary scalar field. More particularly, an exact "normal" 
solution to the static spherically symmetric metagravity equations in the empty, but for a 
singular point, space was written down. The solution is singular in the center and presents 
the black holes filled with graviscalars. It implies the "normal" rotation curves (RC's), 

1. e., those declining asymptotically with distance according to the Newton law. 

In the present paper, a regular "anomalous", missing in GR, solution to the static 
spherically symmetric metagravity equations in empty space is studied. The solution 
naturally results in the "anomalous" , asymptotically fiat RC's. It presents a smooth halo 
as a coherent state of the graviscalar field in the vacuum. Treated in terms of DM the 
halo possesses a finite central density profile reproducing qualitatively the contribution 
to the galaxy RC's due to DM. The way to refine the description of the galaxy halos, as 
composed of the graviscalar DM, is finally indicated. 

2. Anomalous Vacuum Solution 

Unimodular metagravity In the framework of the effective field theory of metric, the 
Lagrangian of the unimodular metagravity looks most generally like:^'^ 



where the graviton and graviscalar Lagrangians Lg and L^., respectively, are as follows: 



L — Lg + Lh + Ljn + Lgh + L 



(1) 




(2) 
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Lh = ^dx-dx-VM, (3) 

with X being the graviscalar fieldS In the above, Kg = l/(87rG)^/2 jg GR mass scale, 
with G standing for the Newton's constant, R is the Ricci scalar, A is the cosmological 
constant, Vh is the graviscalar potential and dx ■ dx = Q^^d^x^uXi with g^j.^ being the 
metric. The Lagrangian for an ordinary matter has some conventional form. In the 
minimal metagravity, we consider, Lg^ = L^h = 0. 

The peculiarity of the graviscalar compared to an ordinary scalar is that the former is 
not independent of the metric, viz., 

X = fl.if. (4) 

^ 9h 

where g = det g^j^i, and gh is a non-dynamical scalar density of the same weight as g. In the 
canonical coordinates, we have gh = —1- The density gt makes x ^ GC scalar and allows 
to bring the theory to the arbitrary observer's coordinates. The parameter Hh stands for 
a unimodular metagravity mass scale additional to the GR Kg. Presumably, Hh < C^(^g)- 
Varying the action S = J d^x^J—gL with respect to g^v, under fixed gn, we arrive at 
the unimodular metagravity equations as follows 

G = Rpiy 'zRgi^iu i^T^u, T^u T\^i/ -\- Tynfjiu Tfi^u, (5) 

with Rfj^y being the Ricci curvature tensor. In the above, T^f^y = Ag^u is the vacuum 
contribution to the total energy- momentum tensor T^^, with Tm^iu being the ordinary 
matter contribution and Thfj,u the graviscalar one. The latter looks like 

Th^,u = d^xduX - Q^X ■ - '^4^5'^^, (6) 

where 

yH = Vh + {dVh/dx + V ■ Vx) , (7) 

with standing for a covariant derivative. The contracted Bianchi identity, V^GJ^ = 0, 
results in the covariant conservation of the total energy- momentum, V^T^^ = 0, instead 
of V ^Tmy = for the ordinary matter alone. (Thereof, the treatment of the graviscalar 
as DM). Eq. ([6]) resembles that for an ordinary scalar field in GR except for the metapo- 
tential Vh superseding the conventional potential Vh. When dealing with the metagravity 
equations, we can proceed in the canonical coordinates, gh = —1, followed by a transfor- 
mation to the observer's coordinates x^. Instead, we proceed directly in x^, with x taken 
as an independent variable, which substitutes a metric element fixed by an additional 
coordinate condition. At that, the unknown gh does not enter the calculations explicitly. 
Having found metric and x we can then through Eq. (jl]) recover in the same coordinates 
the required gh , solving in a sense an inverse problem. 



^The graviscalar being a kind of a "hidden" particle, the related quantities are endowed with a sub- 
script h. 
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In what follows, we restrict ourselves to empty space, Tm/iu = 0. In this case, the 
contracted Bianchi identity results in the graviscalar field equation as follows:^ 

V ■ Vx + dVh/dx = 0, (8) 

with Vfi reduced to 

VH = Vh-Whe-''^^'^. (9) 

Here, Wh is an arbitrary integration constant, not a Lagrangian parameter, which distin- 
guishes the local vacua. 

Particularly, consider the static spherically symmetric configuration of metric and the 
graviscalar field. The line element in the polar coordinates (t, r, 6, (p) looks generally like 

ds'^ = adt^ - bdr^ - cr^dQ, dVt = de^ + sin^ 9d<p^, (10) 

with the three metric potentials a, b and c depending on the radial coordinate r alone. 
The same is supposed about x- Specify r by the coordinate condition ah = 1 and choose 
X as the third independent variable instead of b. Neglect by the potential Vh and the 
cosmological constant A. Putting X = A = a = 1/6 and C = r'^c, we get Eq. as 

{ACXy = ^Ce-^, (11) 

with a prime meaning a derivative with respect to r, and the unimodular metagravity 
equations in the vacuum as^ 



2w 



(CAJ = ^Ce-^, (12) 



3 ^,0 Ku 



(CC")'-^C"' = -^{CX'f, (13) 



2 nl 



(CAJ ~{ACy + 2 = 0. (14) 

By construction, Eqs. f|T2|) - f|T^ are independent, with Eq. ffTTj) being an identity. Instead, 
we choose Eqs. ffTTj) - f|T3|) as the independent ones, with Eq. f|T^ serving as a constraint. 



Anomalous vacuum solution Let first Kh be arbitrary, Kh < 0{ng). At Wh = 0, an 
exact solution to the metagravity equations was given in Ref. 8. At r = 0, the solution 
is singular (refiecting a center point-like matter). This case corresponds to GR in the 
presence of a scalar field. With Wh 7^ 0, there appears a solution regular at r = 0. 
Expanding the unknown functions as the power series in r and equating coefficients at 
equal powers on both sides of equations we have up to terms r^: 
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where ej^ = 2k\/ H?g and = r'^/R'l, with = Gul/wh presenting a characteristic length 
scale squared. Eq. (HM is fulfilled identically up to terms r^. 

Continuing the procedure above we can find the solution with any desired accuracy. 
Namely, knowing the solution in an order r^", n = 0,1,2,... we can first determine X 
from the l.h.s. of Eq. ffTTj) in the next order t'^^"'+^\ Then we can find a and c in the same 
order from Eqs. f|T2l) and f|T3l) . respectively, etc. At that, X{0), a(0) and c(0) are fixed by 
the boundary conditions at r = 0. Both Wh > and Wh < are a priori envisaged. The 
respective solutions are formally related by substitution — r^. For physical reasons, 

Wh> (see later). 

Of special interest is the case ^ 1. Decomposing an exact solution, regular at the 
center, as the power series in Sh (only even powers enter) as X = J^^h^^n, d = J^^'fT^^n 
and c = I]^h"Cn, n = 0,1, . . ., with oq = Cq = 1, we simplify the metagravity equations in 
the appropriate leading orders as follows: 

= Gr'^e-^", (18) 



d 1 


2 dXo 
dr 


d^\ 
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( 2^(11 

I" d7 
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/ ndCl 


dr 



6r^e-^°, (19) 
-l(r'-^)\ (20) 



2 V dr 



with the restriction 



^ - - 2r(ai + Cl)) = 0. (21) 

Clearly, it is possible to add to the solutions for ai and Ci the arbitrary reciprocal terms 
~ 1/r. Assuming no singularity in the center, we omit such contributions. It follows 
from the equations above that the driving term in the system is Xq. Having found the 
latter in a self-consistent manner from Eq. ( fT8|) we can then find Oi and Ci from the two 
other equations with an external source determined by Xq. In particular, it follows that 
ai = Xq modulo a constant which may be put to zero. The leading in Eh part of the 
regular solution given by Eqs. (fTSjl - (fTTI) explicitly satisfies all these equations up to 
accuracy r^. 

To study Xq at > in toto note first of all that there exists an exact exceptional 
solution of Eqs. (USD - ([20]) as follows: 

Xq = ln3r^ (22) 
ai = InSr^, (23) 
Cl = -ln3r2 + 2, (24) 



with the additive constants restricted by the relation ai = Xq and Eq. ( ]2T]) . Present 
further Eq. (fTSi) as follows: 

where Z = Xq — a, with a = InSr^ any real, — oo < cr < +00. Introducing Z = dZ/da 
as an independent variable supplementing Z, reduce the second-order Eq. fl25l) to the 
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equivalent autonomous first-order system: 



da 



Z, (26) 



In the phase plane (Z, Z), there is a single exceptional point Z = Z = 0, defined by the 
requirement dZ/da = dZ /da = 0, other points being normal. Through each normal point 
there should come precisely one phase trajectory {Z{a),Z{a)). The latter ones satisfy 
the equation 

dZ 1,7,1 , , 

with the isoclines dZ/dZ = m being 

(e-^-1), (29) 



2m -\ 

where m is an arbitrary constant. At that, the axes Z = and Z = correspond to 
m —7- ±oo and m = —1/2, respectively. 

Inspection of the phase plane shows that the exceptional point in the center belongs to 
the stable focus type, with all the trajectories winding round the center and approaching 
the latter with a +oo. At that, the exceptional point presents the exceptional solution 
Eq. fl22l) . There is a unique trajectory with Z remaining finite at cr — t- — oo, namely, 
Z — 7- — 1, and behaving thus like Z —a asymptotically. Such a trajectory corresponds 
to the regular at r = solution Xq given by Eq. ( |T5|) with Sh = 0. The rest of trajectories 
satisfy Z — — oo at cr — )■ — oo, with the respective Xq being thus irregular at r = 0. The 
regular anomalous solution is stable against small perturbations of the initial data taken 
on the axis Z = 0, Z < 0, but for very small > 0. In the latter region, the regular 
anomalous solution is to be superseded by a singular normal solution with Vh ^ Rh, where 
Th is the graviscalar radius of a center singularity.^ The account for the latter does not 
significantly affect the halo at r ^ Vh/Rh- 

Altogether, the regular anomalous vacuum solution for the graviscalar field X at ^ 
1 looks like: 

I - Ar^ + + (^(^8)^ at < r < 1, 
°~ \ ln3r2, at r > 1. ^ ^ 

It oscillates around the exceptional solution Xq approaching the latter at r ^ 1. 

3. Anomalous Rotation Curves 

Graviscalar DM The velocity of circular rotation of a test particle in the static spher- 
ically symmetric metric Eq. ( ITOl) is given by 

= TT^- ^31) 

(In r^'c)' 
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So defined velocity transforms as a scalar under the local radial transformations. To get f ^ 
in the leading e/^-order we can put c = Cq = 1. hi Wh> 0, the regular anomalous solution 
results in the scaled RC profile as follows: 

^ elrda, ^ .[r^ - + fr^ + 0(r«), at < r < 1, . . 

'^^^ 2 rfr '^(l, at r > 1, ^ ^ 

where use is made of ai = Xq. At that, the exceptional solution results in the fiat RC 

vlir) = el (33) 

around which all the RC's u^(r), with different Rh, oscillate approaching -u^ at r ^ Rh. 

Let us now interpret RC's in terms of DM. The Newton dynamics in fiat space (a = 
c = 1) with a DM would result in 

where Mh{r) = Att /J^ p/i(r)r^(ir is the DM energy interior to r, with being the DM 
energy density. This implies 

To reproduce the first part of Eq. ( |32l) we should have 

(For a center point-like matter with a — 1~— 1/r and f ^ ~ 1/r, this would give ph = 0.) 
Accounting for Eq. (fTOl) . we get finally the looked- for DM profile as follows (r^ > 0): 

o ~Xo .nJl-r2 + |r4 + C(r6), at < r < 1, 
p,ir) = 2w,e "=^0)^/(3^2)/ at r » 1, ^^^^ 



with the central density 



P.(0) = (3^ 
Rh 



Asymptotically, Mh{r) ~ ef^r/G. Ultimately, such a linear growth should be terminated 
by the potential Vh, which would become significant at the periphery, where Xq gets 
strong. Thus the regular anomalous solution corresponds to a smooth DM halo with the 
finite central density. At that, the exceptional solution results in the cuspy profile 

-P^ir) = (39) 

with the exact Mh{r) = e\r/G. The family of the smooth profiles ph{r), with various Rh, 
oscillates around ph approaching the latter at r ^ i?/^. 

According to Ref. 8, ph = —2Vh may be treated as the energy density of a static 
graviscalar field, incorporating its gravitational energy. This insures a dual field-matter 
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interpretation of the graviscalar halo. In terms of field, the case Wh > presents a lo- 
cal vacuum well, with the metapotential Vh = —Whe~^° due to a coherent state Xq of 
the graviscalar field. In terms of matter, the same case corresponds to the DM distribu- 
tion with ph > which produces in fiat space precisely the same attraction. The case 
Wh < presents a local vacuum bump with the repulsive ph < implying an unstable 
configuration (t>^ < 0). In the GR limit, Wh = 0, the halo clearly disappears. 



Galaxy halos There are numerous studies in astrophysical literature concerning the 
galaxy DM halos. At that, the empirical halo density profiles rely mostly on the two- 
component fits to the galaxy RC's, with the matter and halo contributions added in 
quadrature, f ^ = v"^ + vl, where is a total contribution of the different types of 
luminous matter (disk, gas, bulge) and v"^ is a halo contribution. Thereof, there emerge 
ever growing evidences, based on a vast sample of galaxies of different types, in favour 
of the DM halos with the finite central density profiles (see, e.g., Refs. 9 - 11, with an 
extensive list of references therein). In particular, in Ref. 10 it is found a universal DM 
density profile, extracted from a sample of 36 nearby spiral galaxies, as follows: 

= ,^1^,2 (40) 
1 + [r/RoY 

with po being a central density and a core radius. The empirical smoothness displayed 
by Eq. ( HOl) is at sharp variance with the cuspy form of the cold DM halos. In contrast, the 
vacuum graviscalar halo obtained in the present paper naturally complies with smooth- 
ness. Moreover, Eq. fl37|) closely reproduces the first three terms of the decomposition of 
Eq. pUj) . Nevertheless, there are two differences. First, Eq. (I38l) imphes Ph(0) ~ Rh^^ 
whereas empirically there emerges a constant central surface density of the galaxy halos, 
i.e., PqRq ~ const. Second, Eq. f p7|) is three times lower asymptotically compared to 
Eq. di]). 

The reason of the discrepancy may be as follows. We have restricted ourselves by the 
simplest model with the spherical graviscalar halo in empty space. Such a halo may serve 
just as a prototype for the real galaxies. To confront the theory with the data the luminous 
matter should also be accounted for. This would give v"^ = vf^ + v^^, with v^^ being the 
effective graviscalar halo contribution in the presence of matter. Because of a coherent 
nature of halo the deformation of the latter, both in magnitude and sphericity, may be 
significant within the region of intersection of matter and halo. It is rather f^^, to which 
the two-component fit Eq. fHUl) is to be applied, than due the vacuum graviscalar halo. 
To distil the galaxy RC sample from the matter contribution as far as possible, live aside 
the points from Ref. 10 which correspond explicitly to the luminous matter dominance at 
the distances at hand. Inspection shows that there are at least four, out of 36, points on 
the log Po - log Rq plot, with the extremely large Ro, to be dropped off. (Incidentally, such 
pecuhar Rq are strongly model dependent.) The rest of points lies much more compactly, 
the subsequent results being less sensitive to further reducing the galaxy sample. On 
the reduced sample, the dependence logpo ~ ~ log-Ro looks less prominent. On the other 
hand, a much wider sample of galaxies of different types still supports the constant central 
surface density law.^^0 To settle the question in the metagravity framework the account 



"The latter data may also indicate some deviations from Eq. pO|) both at very small and very large 
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for the luminous matter, with an axisymmetric distribution resulting in the graviscalar 
halo asphericity, is required. 

Finally note that it is the Lagrangian parameter Sh = \/2Kh/Kg, which sets the scale 
of the asymptotically flat RC's due to the graviscalar halo in the vacuum. So, taking for 
galaxies asymptotically Vh{oo) ~ 100 km/s we would expect that Eh = fh(oo)/c ~ 10^^. 
With Hg = 2.4 X 10^^ GeV, the unimodular metagravity mass scale, Hh ~ 10^^ GeV, would 
approach the GUT mass scale. 



4. Conclusion 

The regular anomalous solution to the static spherically symmetric metagravity equations 
in empty space presents a viable prototype model for the smooth galaxy halos character- 
ized by the finite central density profiles. It goes without saying that once the solution 
models the descent density profiles of the DM halos, it provides to the same extent all 
the other gravitational effects of such halos. The hypothesis about the graviscalar origin 
of DM in the framework of the unimodular metagravity finds thus its preliminary confir- 
mation. To further verify the theory studying the graviscalar halos in the presence of the 
axisymmetric matter distribution is in order. 
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